MIT-CTP 3881 



A calculation of the bulk viscosity in SU(3) gluodynamics 
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We perform a lattice Monte-Carlo calculation of the trace-anomaly two-point function at finite 
temperature in the SU(3) gauge theory. We obtain the long-distance properties of the correlator 
in the continuum limit and extract the bulk viscosity £ via a Kubo formula. Unlike the tensor 
correlator relevant to the shear viscosity, the scalar correlator depends strongly on temperature. If 
s is the entropy density, we find that f/s becomes rapidly small at high T, C/s < 0.15 at l.65T c and 
C/s < 0.015 at 3.2T C . However C/s rises dramatically just above T c , with 0.5 < C/s < 2.0 at 1.02T C . 

PACS numbers: 12.38.Gc, 12.38.Mh, 25.75.-q 



Given the success of the hydrodynamical descrip- 
tion [l[ of high- energy heavy ion reactions, it is of primary 
interest to compute the shear and bulk viscosities of the 
quark-gluon plasma in the region T c < T < 3T C relevant 
to the RHIC and the forthcoming LHC experiments. 

Phenomenology combined with hydrodynamics calcu- 
lations point to a very small shear viscosity to entropy 
density ratio, rj/s < 0.2 _m, and viscous calculations 
find even smaller values 3j. The main uncertainty in 
these determinations is the sensitivity to the initial con- 
ditions of the hydro-regime. A recent lattice calculation 
of -q/s in the purely gluonic plasma [1] found values in 
the range 0.1-0.2, under an assumption that we shall 
return to below. These values are significantly smaller 
than the leading-order perturbative result r)/s f=a 0.8 for 
a s = 0.3|. 

In this Letter we carry out a lattice calculation of the 
bulk viscosity, £, in the SU(3) gauge theory, reusing much 
of the technology developed in Q . The calculation relies 
on a Kubo formula and on the reconstruction of the real- 
time spectral function 0, 0, HI from the Euclidean trace- 
anomaly correlator. This is known to be a numerically 
hard problem, in particular some assumptions have to 
be made on the spectral function in order to extract the 
viscosity — typically, an assumption of smoothness. In 
this respect, the bulk viscosity is a more favorable case 
than the shear viscosity at high temperatures, because 
the spectral function is rigorously known to be smooth at 
T — > oo. This is related to the fact that the bulk viscosity 
vanishes for any conformal field theory, and non-Abelian 
gauge theories become conformal at high energies by be- 
coming free. 

It has recently been argued [9| , based on an exact sum 
rule, lattice data on e — 3P and a simple model for the 
spectral function, that the bulk viscosity rises sharply 
close to the (first order) phase transition. This is remar- 
kable, because a gas of structureless point particles has 
negligible bulk viscosity both in the non-relativistic and 



the extreme-relativistic limits [Io| . and the bulk visco- 
sity is often neglected for that reason even when the shear 
viscosity is taken into account. Furthermore, the leading- 
order perturbative QCD result |ll| for C/s is numerically 
very small: to within 10% for 0.06 < a s < 0.3, 



C/s = 0.020a 2 s (N{ = 0). 
It satisfies the relation [l(| C ~ l^wih -~ v 1] 



(1) 

, being the 

velocity of sound in the medium. Instead, AdS/CFT cal- 
culations that incorporate deviations from conformality 
at strong coupling find C oc rj[^— v^] [l2j. These formulas 
suggest that the dynamics responsible for the bulk visco- 
sity away from conformality, v 2 s < | , is different at strong 
and weak coupling. 

We will show that C/s depends strongly on tempera- 
ture: it is phenomenologically negligible above 2T C , while 
it becomes 0(1) just above T c . We will also acquire some 
knowledge of the spectral function from first principles 
thay may help interpret this striking behavior. 



Methodology 

The explicit form of the Euclidean ener gy- m omentum 
tensor can be found for instance in [181 ]: we write 
Too = T'oo + \6, where 9 is the trace anomaly. With 
L = l/T the inverse temperature, we consider the two- 
point functions (0 < xq < Lq) 

C e (x Q ) = L 5 fd 3 x (0(0)0(^0, x)) (2) 
C- (x ) = L 5 fd 3 x (T o(0)Too(£o,x)}, (3) 
C'oe(xo) = L 5 Jd 3 x (T O o(O)0(xo,x)). (4) 
This correlator is represented by a spectral function, 
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in terms of which the bulk viscosity is given by [if 



'Electronic address: meyerh@mit.edu 



(5) 



(0) 



j=0 



2 



C(x ) 



10 - 



1.02T C ,LT=3 
1.24T c ,LT=5/2 
1.65T C , LT=5 
2.22T C , LT=6 
3.22T, LT=6 



T x n 



0.2 



0.3 



0.4 



0.5 




0.005 0.01 



0.015 



0.02 



0.025 



FIG. 1: The correlator Cb(xq) on Lo/a — 8 lattices. 



Important properties of pg are its positivity, pg{us)/uj > 
and parity, pg(-uj) = —pg{u). At tree level, Cg (xq) — 
^ a 2 s N 2 [f(r) - fl], with r, d A and /(r) as in 0, and 



tanh jujLq ' 



A 



t.i. 



lla s iV 
3(4tt) 2 



(7) 



By contrast with the shear viscosity case [3] , the absence 
of a 5-function at the origin corresponds to the fact that 
the theory is conformal at tree-level and hence the bulk 
viscosity vanishes. 

Because / d 3 x(T 00 (x)O(0)} = T 2 d T {0) T for any local 
operator O and xo ^ 0, the two other correlators are 
related to Cg by 



be(T) 



±Cg(x ) + b (T), 
-\Cg{x Q )+bg{Tl 



l T - 
2^ 



3 <9 T (e-3P) 
T- 3 d T (e-3P). 



3T~ 



(8) 
(9) 
(10) 

(11) 



i.e. they only differ by an cco-independent shift. In per- 
turbation theory, b 9 (T) = ^^ s ) 2 + 0(a 3 ) and b„(T) 



is also easily extracted from [17j. This completes the 
perturbative result for the Euclidean correlators up to 
O(al). 

We finally introduce the following moments of the spec- 
tral function (n = 0, 1, . . .): 



2)i 



d 2n Cg 



sinh ujLq/2 



dxl n 



(12) 



x Q =L /2 



The latter equality implies that they are directly acces- 
sible to lattice calculations. For a spectral function that 
grows like lo a at large frequencies, it is natural to quote 



the pure numbers 
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FIG. 2: Continuum extrapolation of the first two moments 
(uj°) and (uj 2 )/(12QT 2 ) of pe{oj,T) for T = 1.24 and 1.65T C . 
The (to ) continuum limits are resp. 16.3(5) and 8.5(8). 



Numerical results 

We simulate the isotropic Wilson action with a two- 
level algorithm [14|, as in [J], and emplo y t he 'bare- 
plaquette' discretization of Too and 9 as in [18|. We use 
the data of [l3| to determine the value of Tro and to com- 
pute the lattice beta-function, and use T c r = 0.75(1) [HI 
to obtain T/T c . 

We present the Cg correlator, computed on L /a = 8 
lattices in the temperature range T c < T < 3T C , in Fig.[T] 
It has been 'improved' by the same technique as in [J], 
so that cutoff effects are 0(g 2 a 2 ). The deviations from 
conformality are significantly larger than those seen in 
the tensor channel relevant to the shear viscosity As 
T approaches T c from above, they become very large. 

We investigated finite- volume effects at 1.65 and 1.24T C 
by varying LT from 3.5 to 5 and 2.5 to 3.67 respectively; 
no statistically significant dependence was observed. 

Figure [5] shows the continuum extrapolation of the first 
two moments of pg from lattices Lo/a = 6, 8, 10, 12. The 
scaling region for starts at Lo/a = 6, while for (u> 2 ) 
it presumably only starts at Lo/a = 8 or 10. We also 
find s/T 3 = 4.46(6) and 5.48(7) at T = 1.24 and 1.65T C 
in the continuum. 

On the lattice, the relations (P fTTj) are only satisfied up 
to 0(a 2 ) corrections. They can be used to over-constrain 
the continuum extrapolation of the Cg correlator. It 
turns out that the direct measurement of Cg is more ac- 
curate than the estimators based on relations (|8l9p by at 
least a factor five, and that Cg has large discretization er- 
rors; we have nonetheless checked that the data is consis- 
tent with relations (f51 TTTj) in the continuum limit. While 
the spectral representation ([5]) is in principle only valid 
in that limit, smaller lattice spacings will be necessary to 
extrapolate the whole correlator to the continuum. 

We thus proceed to reconstruct the spectral function 
at the smallest lattice spacing, based on N points of 
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FIG. 3: The result for p e {uj) from L /a = 12 lattices. The 
bulk viscosity is given by C,/T A — (7r/18) x intercept. The 
oscillating curve is the (rescaled) resolution function 5(0, ui). 

the correlator at abscissas {xq > 2a}, with (strictly 
positive-definite) covariance matrix Sw . We write [4] 
p(u>) — m(u>) [1 + a(u)) where m(uj) > has the high- 
frequency behavior of Eq. [7J In practice we have used 
m(uj) = A t .i.w 4 /[tanh(io;Lo) tanh 2 (co;Lo)] where c is a 
tunable parameter typically set to 4. We build an esti- 
mator of a(ui), a(uj) = ^2^ = iC£Ui(uj), using the same ba- 
sis {ug} as in [4|. Upon discretization of the u> variable, 
M(xq,lo) = K(xQ,oS)m((jS) becomes an N x N u matrix 
Mij, and the functions {ue}f =1 are the columns of the 
matrix U in the singular-value decomposition of M t = 
UwV t . The positivity of p(oj) = m(u))[l + a(oj)] pro- 
vides an important a posteriori consistency check. The 
spectral functions obtained in this way from Lq/cl = 12 
lattices are shown on Fig. O for which N = 4. 

Solving for the coefficients ce results in a(uf) = 
E^=i A Ci«H, with AC, = C(4 4) ) - jduMiiw) and 
qi(uj) — Vi^wJ ui(ui), which is related to the genuine 
spectral function by 

a(cu) = f™du' 6(uj,w') a(u'), (13) 
6{u,u') = Zti 9iM Mi(cj'). (14) 

The resolution function 5(0, to) is plotted in Fig. [3] 
Clearly it has the effect of smoothing a(u>) over a region 
Aui « 5T around the origin; note that a(u>) "^>°° o. 

In solving inverse problems, it is often found neces- 
sary, when N becomes large, to introduce a 'regulator', 
whose role is to prevent the reconstruction from becom- 
ing unstable, at the cost of making the resolution function 
S(uj, ui') less sharply peaked around ui' = u. One way (2p| 
to implement this in our method is to replace the qi(ui) 
by [wV t + \V t S]~ 1 a ug(ijj), where A is typically chosen 
so that p yields a x 2 /d.o.f. of order unity. We have how- 
ever not found it necessary to make A ^ in any of the 
cases presented here. 



FIG. 4: The bulk viscosity in units of the entropy density, as 
a function of the conformality measure e ~^f . The statistical 
errors are shown, as well as the bounds explained in the text. 
The solid line is the perturbative prediction Eq. [T] naively 
continued beyon a 3 — 0.3 by the dashed line. 

The value of p/uj\ui=o and its statistical error, together 
with the resolution function, form an assumption-free re- 
presentation of the information on ( contained in the 
lattice data. In order to quote a systematic error on £, 
some modelling of the spectral function is necessary. 

Take for instance T — 1.24T C , where our best estimate 
is £/s = 0.065(17) s tat- Increasing c from j to | reduces 

C/s to 0.010(2) but causes p to assume negative values 
around u> — 10T by one standard deviation. One may 
therefore regard this as a lower bound on £/s. We also 
infer a conservative upper bound based on the positivity 
of p 0| • Setting p = p t .\. for ui > 20T C , and introducing 
a Breit-Wigner of width 2T at the origin, we find that 
C/s < 0.37 at 90% confidence level in order for that con- 
tribution to the correlator not to exceed Cq(Lq/2). We 
obtain similarly C/s < 0.15 at 1.65T C . 



Bulk viscosity & non-conformality 

It is clear from Fig. [3] that for T ^ T c the spectral 
function increases overall and also the spectral weight 
clusters around the origin. Both features contribute to 
an increase in the bulk viscosity, and they are robust 
predictions, since they essentially correspond to the in- 
crease of (w°) and the decrease of (w 2 )/(T 2 {w )). We 
now study the relation of this effect to non-conformality 
in more detail. As a measure of (non)-conformality, we 
use e ~+p , which is easier to compute than the velocity of 
sound but to which it is related at weak coupling by 

1 - 3v 2 s = l(e - 3P)/(e + P) [1 + 0(a.)]. 

We obtain e — 3P from [16| by quadratic interpolation 
and compute e + P ourselves. In Fig. [H we show the 
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ratio C/s obtained on L^/a — 8 lattices, where we have 
data at five different temperatures, as a function of this 
measure. We define a renormalized, "(e — 3P)-scheme 



coupling", via 



e-3P 
c+P 



55N^a 2 s (27TT) 

967T 2 



and insert it into 



Eq.[TJ That leads to the curve shown on the same plot. It 
lies generally below the lattice points and, unsurprisingly, 
does not account for the sharp rise of C close to T c . 



Conclusion 

We have computed the correlation functions of the 
energy-momentum tensor relevant to the bulk viscosity 
in the SU(3) pure gauge theory, as well as the entropy 
density s, to high accuracy. Our best estimates of the 
bulk viscosity, obtained on Lq/q = 12 lattices, are 



0.008(7) 
0.065(17) 



0.15 
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(T = 1.65T C , LT = 
(T = 1.24T C , LT = |). 



(15) 

where the statistical error is given and the square bracket 
specifies conservative upper and lower bounds. This con- 
stitutes the first successful lattice calculation of the bulk 
viscosity of an SU(AT) gauge theory. In Eq. [15j the 
bounds hold in the continuum. Since we now have esti- 
mates of both r\ [i| and ( at 1.24T C , we can ask whether 
the relation ( w 15f?[± -v 2 ] 2 holds. Using v 2 s « 0.25 [HI, 
the formula predicts (,jr\ w 0.1, while our central value is 
0.6. Our data thus points to a much larger bulk viscosity 
than the formula suggests. A non-negligible bulk visco- 
sity is a signature of the interplay between translational 
and internal degrees of freedom; a small shear viscosity 



that of a strongly interacting system. 

The bulk viscosity is non- vanishing below 1.3T C , where 
deviations from conformality in e and P are large. In 
particular we observe a sharp rise in £ in an Lo/a = 8 
simulation just above the first order phase transition, 



C/s = 0.73(3) 



2.0 
0.5 



(T = 1.02T C , LT = 3). (16) 



Looking to 'ordinary' substances, it is worth noting that 
the bulk viscosity of nitrogen remains finite on the liquid- 
vapor coexistence line [2l|, while it becomes very large 
close to the endpoint (for the case of xenon see [22|). Al- 
though QCD with physical u, d, s quark masses seems to 
undergo a crossover rather than a phase transition [l9l |. 
the rise of £ is related to the large rate of increase of 
(e — 3P)/T 4 [H and to the clustering of the spectral 
weight towards u> — (Fig. [3]). We therefore also ex- 
pect an increase in the bulk viscosity around 200MeV in 
the real world. The effect could manifest itself in heavy 
ion collisions by an enhanced entropy production in the 
later stages of the reaction, provided the viscous hydro- 
dynamical description is still valid at that time. 
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